Revised Iterative Solution of 
Ground State of Double- Well Potential 



Zhao Wei-Qin^' ^ 



1. China Center of Advanced Science and Technology (CCAST) 
(World Lab.), P.O. Box 8730, Beijing 100080, China 
2. Institute of High Energy Physics, Chinese Academy of Sciences, 
P. O. Box 918(4-1), Beijing 100039, China 



Abstract 

A revised new iterative method based on Green function defined by quadratures 
along a single trajectory is developed and applied to solve the ground state of the 
double-well potential. The result is compared to the one based on the original 
iterative method. The limitation of the asymptotic expansion is also discussed. 
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1. Introduction 



The double-well potential in one dimension, 

V{x) = \g\x^ - l)^ (1.1) 

is specially interesting since it has degenerate minima and could be served as a simple 
example of the bound-state tunnelling problem in quantum mechanics. However, it is 
a non-perturbative problem to solve the Schrocdinger equation for this potential due to 
the tunnelling effect between the two minima. The asymptotic series of the ground state 
energy is meaningful only for quite large g (e.g. fif > 6). Even for such large g the obtained 
plateau in the energy expansion series does not necessarily consistent to the exact solution. 
An effective method to obtain the convergent solution of this problem for any values of g 
is needed. 

Recently an iterative solution of the ground state for the one dimensional double- well 
potential is obtained[l] based on the Green function method developed in ref.[2]. This 
Green function is defined along a single trajectory, from which the ground state wave 
function in N-dimension can be expressed by quadratures along the single trajectory. 
This makes it possible to develop an iterative method to obtain the ground state wave 
function, starting from a properly chosen trial function. The convergence of the iterative 
solution very much depends on the choice of the trial function[l]. 

However, in the original iterative solution the obtained correction for the trial wave 
function is in the form of a power expansion, while the bound state solution should be 
in the form of an exponential when the coordinate variable approaches infinity. Recently 
a new revised iterative procedure [3] based on the same Green function is developed and 
applied to solve anhamornic oscillator and Stark effect. This method has some advantages 
compared to the original one. It not only gives an exponential form for the correction of 
the trial wave function, but also spends much less time in each iteration due to less folds 
of integration. It is natural to try this method to solve the double-well potential which 
has been attracted much attention. 

In Section 2, a brief introduction is given about the Green function method based on 
the single trajectory quadrature. Special discussion is given to the revision of the iterative 
formula. The revised iterative formula for the double-well potential is given in Section 3, 
together with the trial function and the boundary condition for the lowest even state. The 
numerical results based on the revised iterative formula and the original iterative method 
are collected in Section 4, together with a comparison to the asymptotic expansion. A 
compact expression of the asymptotic expansion is derived in Appendix A. Finally some 
discussions are given at the end. 
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2. Green Function and the Revised Iterative Solution 



For a particle with unit mass, moving in an N-dimensional unperturbed potential Vo(q), 
the ground state wave function $(q) satisfies the following Schroedinger equation: 

i/$(q) = E$(q), 

where 

H = r + yo(q) = -lv' + yo(q). 

Assume the solution of eq.(2.1) could be expressed as 

$(q) = e-^('i). 
Introduce a perturbed potential U (q) and assume 

^q) = W + t^(q) 

n = r + y(q)^_lv2 + F(q). 
Define another wave function ^'(q) satisfying the Schroedinger equation 

£ = £; + A. 

Let 

The equation for r and A could be derived easily [2]: 

VS-VT + \[{VTf - VV] = {U- A). 
Consider the coordinate transformation 

?2, ^3, • • • , (In) {S, ai, 02, • • • , aN-i) = {S, a) (2.10) 

with a — (cKi, a^, - • ■ , ckat-i) denoting the set oi N — 1 orthogonal angular coordinates 
satisfying the condition 

VS-Vai = 0, (2.11) 
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for i = 1, 2, • • • , — 1. Some useful quantities for this coordinate transformation [2] are 
listed in Appendix B. Similarly to the discussions in Ref.[2], introducing the ^-function 
in S-space: 

and define 

C = e[{VSy]-^ = ehl. (2.13) 

Using 

V5-VC=1, (2.14) 
it is easy to derive the following equation 

r = (1 + CT)-'C[{U - A) - \{VTf]. (2.15) 

When the N-dimensional variable q is transformed into {S, a), T — — could be 
decomposed into two parts: 

T^Ts + T^, (2.16) 

where Ts and consist only the differentiation to S and a, respectively. The detailed 
expression of Ts and Tq, is given in Appendix B. Now another Green function could be 
defined as [2] 

D = -2ee^'^ee-^'hsK (2.17) 

and it is related to C in the following way [2] : 

(1 + DT^Y^D = (1 + CT)-^C. (2.18) 
Therefore, from (2.15), we have 

T = (1 + DT^)-'D[{U - A) - l(Vr)2]. (2.19) 
The exphcit expression of r based on (2.19) is 

r = -2 [\^s'^dS' r e-^^"hs"hadS"{l+T^D)-'[{U-A)-hvrf]. (2.20) 
Jo ha Jo 2 

Therefore, we have 

= r e-'''hs'KdS'{l+TaD)-\U - A) - Uvrf]. (2.21) 
Jo I 



1 ha _25 5t(5', a) 



2hs dS 



The left hand side of Eq.(2.20) approaches to when 5" — > oo, so is the right hand side, 
i.e., 

/ e-^%shadS{l + T^D)-'[{U - A) - -{Vrf] = 0, (2.22) 
Jo I 

which is correct for all a. Integrating over da. — ^J^dai and because of 

hshaTa = X] j da^shaTaT = 0, (2.23) 

j=i j j j 

we derive 

/' hsh^dadSe-^'^liU - A) - -{Vrf] = 0. (2.24) 

Denoting dq = hghadadS , from Eq.(2.3), we reach a new expression of the perturbative 
energy 

Idc^^[U-UVrn 

Based on Eqs. (2.25) and (2.15) or (2.19) we have the new iteration series 
A„, = 



/dq$2[C/-i(Vr„_i)^] 



Tn = {l + DT^y'D[{U-A^)-^{VTn-if] (2.26) 
or 

Tn = (l + CT)-iC[(t/-A0-^(Vr„_i)2]. 

For later convenience this new iteration series is name as the r-iteration in this paper. 
Let us compare the above r-iteration with the original one derived from the equation for 
/ = e^'^ and A in Ref.[2], which is named as /-iteration in this paper: 

^ Jdq ^'U fn-l 

fn = l + {l + DT^)-'D{-U + An)fn-i (2.27) 
or 

fn = 1 + {1 + CT)-'C{-U + An)fn-1. 

There are several advantages for the r-iteration: 

1) It directly gives an exponential form for the perturbed wave function e""^. This 
result is consistent with those obtained using the series expansion of {Si} and {Ej} (See 
Section 1 of Rcf.[l]); 

2) The iteration process in this formula is more transparent; 

3) The calculation of the perturbation energy is much simpler. 
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3. Revised Iterative Formula for the Double-well 
Potential 



In this section the revised iterative formula is apphed to solve the ground state for the 
double-well potential. For the Hamiltonian H = T + V let us introduce the wave function 
for the lowest even eigenstate as ipgy, satisfying 

Hlpe^ = Ee^lpe^. (3.1) 

In the following we are going to introduce the trial wave function (p^^ for this state, 
satisfying 

{H + Wev)(pev = {Eev + £ev)<Pev = g<Pev (3.2) 

Assuming 

ll^ev = 0e^;e"'^'^^ (3.3) 

the final energy and wave function {ipev, E^v} could be obtained by solving the correspond- 
ing equation of {rev, ^ev} based on the revised iteration method introduced in Section 2. 
The key to choose the proper trial function is to satisfy the necessary boundary conditions 
of the state. For the ground state, namely the lowest even state, we have 

'4'ev{-x) = Ipevix) 

i^'eM = (3.4) 
ilev{oo) = 0. 

The trial wave function should satisfy similar boundary conditions, namely. 



■'ev 



-X) = (j)ev{x) 
0e.(O) = (3.5) 

>.Joo) = 0. 
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Following the steps in Section 2 of Ref.[l], we introduce, for x >0, 

gSoix) = |(x-lf(x + 2), 

S,ix) = In^, (3.6) 

(p+ix) = (f>+{-x) = e-^^°(^)-^i(^) = e-^^°(^)(^— ), 

1 I X 

o 

0_(a;) = e~^'^°^~^^~'^^^^^ = e~3^e"'"^'^°*-'^'*(— — ). 

It is easy to see that satisfies 

{T + V + u)<P+ = g<P+, (3.7) 

where 

nx) = ^/(x^-1)^ 

(1 + xy 

Following the necessary boundary conditions (3.5), we could choose the trial function as 

, , i (p+ix) + ^(p^ix), for 0<x<l ^ ^ 

0e.(x) = 0e.(-x) ^ I (1 + lie-h)M^), for x > 1. ('-^^ 

It is easy to proof that the above trial function satisfies Eq.(3.2) with 

Wev{x) = Wev{—x)=u{x)+gev{x), for x>0 (3.10) 

... 2g to-l)e----°;-^Y- ^ forO<X<l 

for X > 1. 

Although the potentials V + Wev is not continuous at a; = 1, it could be proved that (pev 
and are continuous at any x, including x — 1 and x — 0. 
Now introduce the ^-function in x-space: 

and define the Green function 

D,, = -29<P:,W,,. (3.12) 
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Following similar steps as in Section 2, we obtain the equations for {tsv, Eev} as follows: 

Tev = De^[-{Wev- ^ev) -]j^{'r'^v)\ (3-13) 



oo 




oo 




The solution of the ground state is 

i^ev = (Peve~''^"', Eev^g-Sev (3-14) 

Similar procedure could also be taken by introducing 

D+ = -264)^9(1)% (3.15) 

and 

r+ = D4-{u-8^)-\{T'^n (3.16) 

oo 

p - 

! <t>+ dx 



It should be noticed that although if)^ = 0+ e~'^+ does satisfy the equation 

= (3.17) 

with = g — the solution ■0+ is not the eigenstate of the Hamiltonian H — T + V, 
because the trial function 0+, as well as the solution ip+ do not satisfy the necessary 
boundary conditions. However, we still like to keep this solution here since an analytic 
expression of and ip^ has been obtained in terms of the asymptotic power expansion 
of 1/gf in Ref.[l]. This makes it possible to check the accuracy of the iteration procedure. 
In Appendix A compact expressions of the asymptotic power expansion of and are 
derived. 

Based on the explicit integral expression of D^^ and taking derivatives of Tgy and 
T+ in the first equations of (3.13) and (3.16), we obtain 

rM = 20-'(^) J<PlM[My) - U + l^Mf] dy, (3.18) 



oo 

j<i>i{x) (wux) + UK.{x)r)dx 



£ 







oo 

/ ^IvW) dx 
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and 



^[{x) = 20-^(x) / <l>Uy)My) - + 2«(y))1 dy, (3.19) 



oo 

/ (\)\{x) dx 



The two sets of equations (3.18) and (3.19) are for the pairs of Tg^, £ev and t+, £+, which 
could be solved iteratively in the following way: Introducing the initial conditions 

4.,o = 0, r^,,o = (3.20) 

and 

£+,0 = 0, r;o = 0, (3.21) 

we have, for n = 1, 



CXD 



/ (ply{x) Wev{x)dx 
^ev,i — 55 (3.22) 





X 





oo 



and 



/ 0+(3^) u{x)dx 

£:+,! = . (3.23) 

! (p\{x)dx 







t;i(x) = 2r+\^) j 4>l{y){u{y) - £^,,)dy 



oo 

= -20;2(^) j ct>l{y){u{y) - £+,^)dy 



For n > 1 we have 



9 



£ev,n — £ev,l H 55 (3.24) 





X 



and 



= W(^) + 2</'e^^(^) / <^L(y)[(^eM-<^e^,n) + ^(Te^,n-l(y))Vz/ 



oo ^ 

= T^,,i(x) - 2<t>:^{x) I - «^e.,n) + ^{rL,n-i{y)y]dy 



£+,1 + ^ ^ . (3.25) 

/ (t)\{x) dx 





a; 



<,n(^) = <,i(^)+20;^(^)/0l(y)[(^+,i-^+,n) + ^(rVi(y))Vy 



oo 

X 

For comparison we list in the following also the /-iteration formula for the double- well 
potential based on (2.27) [2]: 

oo oo 

fn(x) = 1 - 2 / r\y)dy 1 0'(^)(«;(^) - En)fn-i{z)dz (3.26) 

X y 

oo oo 

En ^ (f)'^{x)fn-l{x)w{x)dx / / (Z)"^ {x) fn-l{x)dx, 



where and could be either for the even ground state "ev" or the "+" state. 

In the following section we are going to show the iterative results of and Sev,n 
based on the r-iteration formula Eqs.(3.24) and (3.25), together with the result from 
the /-iterative formula (3.26). The obtained „ is also compared to the asymptotic 
expansion for different g, obtained from the compact expression in Appendix A. 
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4. Numerical Results and Discussions 



Now let us look at the results based on the r-iteration formula (3.24) and (3.25). Starting 
from the trial wave functions (pev and (j)^ defined in (3.9) and (3.6), the energies E^v and 
after the first 4 steps of iteration with 5-fold integrations are listed in Tables 1 and 2, 
together with the energies i?et),5 = 9 ^ ^eu.s and i^+^s = 9 — ^+,5- In Ref.[l] a convergent 
iteration method to solve the lowest states of the double-well potential, i.e. the /-iteration, 
has been given. However, no numerical results have been provided. For comparison the 
numerical results for S^v and £+ based on the /-iteration method is also listed in Tables 1 
and 2. For the same folds of integrations the solution of the /-iteration could only reach a 
lower accuracy, comparing to the r-iteration, because each step iteration of r' needs only 
one fold of integration, while each step iteration of / depends on 2-fold integrations. The 
accuracy of the obtained energies is higher when g becomes larger. 



Table 1. £ev,n and Eev^ = 9 - ^ev,5 



9 


n 


1 


2 


3 


4 


5 




0.05 


r-iter. 


-0.0341 


-0.0172 


-0.0158 


-0.0163 


-0.0164 


0.0664 


0.1 


r-iter. 


-0.0118 


-0.0022 


-0.0016 


-0.0017 


-0.0017 


0.1017 


0.3 


r-iter. 


0.0963 


0.0973 


0.0973 


0.0973 


0.0973 


0.2027 


0.5 


r-iter. 


0.2035 


0.2060 


0.2060 


0.2060 


0.2060 


0.2940 


0.5 


/-iter. 


0.2036 


0.2056 


0.2060 








1 


r-itcr. 


0.4135 


0.4310 


0.4312 


0.4311 


0.4311 


0.5689 


1 


/-iter. 


0.4135 


0.4267 


0.4302 








3 


r-iter. 


0.4757 


0.5105 


0.5166 


0.5173 


0.5173 


2.4827 


3 


/-iter. 


0.4757 


0.5053 


0.5141 








6 


r-iter. 


0.29204 


0.29399 


0.29420 


0.29422 


0.29422 


5.70578 


6 


/-iter. 


0.29204 


0.29393 


0.29419 








7 


r-itcr. 


0.27884 


0.27957 


0.27962 


0.27963 


0.27963 


6.72037 


7 


/-iter. 


0.27884 


0.27955 


0.27961 








8 


r-iter. 


0.27231 


0.27265 


0.27266 


0.27266 


0.27266 


7.72734 


8 


/-iter. 


0.27231 


0.27261 


0.27266 
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Table 2. S+^n and .g+^s - g+,5 



9 


n 


1 


2 


3 


4 


5 


^+,5 


1 


r-iter. 


0.4135 


0.4310 


0.4312 


0.4311 


0.4311 


0.5689 


1 


/-iter. 


0.4135 


0.4267 


0.4302 








3 


r-iter. 


0.3221 


0.3257 


0.3258 


0.3258 


0.3258 


2.6742 


3 


/-iter. 


0.3221 


0.3254 


0.3257 








6 


r-iter. 


0.27989 


0.28040 


0.28041 


0.28041 


0.28041 


5.71959 


6 


/-iter. 


0.27989 


0.28039 


0.28041 








7 


r-iter. 


0.27461 


0.27494 


0.27494 


0.27494 


0.27494 


6.72506 


7 


/-iter. 


0.27461 


0.27493 


0.27494 








8 


r-iter. 


0.27091 


0.27113 


0.27113 


0.27113 


0.27113 


7.72887 


8 


/-iter. 


0.27091 


0.27113 


0.27113 









It is shown clearly that the obtained energies E^.^ and Ej^ are lower than g and E^.^ < 
which are reasonable. When g increases the two energies become very close to each 
other and the second step of iteration gives already quite accurate result. It is interesting 
to notice that for g = 1 the trial function 0+ and (pev are the same, therefore starting 
from the obtained -0+ by iteration is the exact ground state wave function for gf = 1. 



Table 3. £ev and Eev — g — £( 



9 


'-'ev 


E 

-'-'ev 


9 


'-'ev 


E 

-'-'ev 


0.05 


-0.0164 


0.0664 


2.2 


0.5951 


1.6049 


0.1 


-0.0017 


0.1017 


2.5 


0.5745 


1.9255 


0.3 


0.0973 


0.2027 


2.7 


0.5539 


2.1461 


0.5 


0.2060 


0.2940 


3.0 


0.5173 


2.4827 


0.7 


0.3065 


0.3935 


4.0 


0.3984 


3.6016 


1.0 


0.4311 


0.5689 


5.0 


0.3273 


4.6727 


1.5 


0.5598 


0.9402 


6.0 


0.29422 


5.70578 


1.7 


0.5851 


1.1149 


7.0 


0.27963 


6.72037 


2.0 


0.5990 


1.4010 


8.0 


0.27266 


7.72734 



In Table 3 listed are the obtained Sev for different g based on the r-iteration. The 
obtained Sev is negative for very small g and increases when g increases from 0.05 to 
2.0, then decreases when g increases further. However, the energy for the ground state 
Eev — g ~ ^ev monotonically increases with increasing g. 
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N 

Table 4. 8+, = g — Sj^ and = J2 ^m+i/g 

rn=0 



9 




E+ 






El 


0(e-t5) 


1 


0.4311 


0.5689 


- 


- 


- 


~ 0.3 


2 


0.3664 


1.6336 








~ 0.07 


3 


0.3258 


2.6742 








~ 0.02 


4 


0.3024 


3.6976 








~ 0.005 


5 


0.2888 


4.7112 








~ 0.001 


6 


0.28041 


5.71959 


10 


17 


0.2807 


~ 3 X lO-'* 


7 


0.27494 


6.72506 


14 


22 


0.27501 


~ 4 X 10-5 


8 


0.27113 


7.72887 


12 


34 


0.27115 


~ 2 X 10-5 


9 


0.268336 


8.731664 


11 


35 


0.268339 


~ 6 X 10-6 



The energy and are given in Table 4 together with e+ calculated from the 
asymptotic power expansion to Xjg. It can be seen that only when g is large enough, 
say 9' > 6, the asymptotic expansion is meaningful. For such large g the obtained 
from the iteration and from the power expansion are comparable up to a quite accurate 
level. The power expansion oi \,o Xj g is an asymptotic one. For a fixed and large 
enough (7-value the summation up to a certain number of terms becomes stable. When 
increasing the number of summed terms further a plateau of the energy E^ (within the 
number of summed terms A^rnin < N < N^an shown in Table 4) is obtained, which gives 
the i?+-value accurate to a certain level. However, beyond a certain number of terms 
{N > Nmax) the result becomes unstable and meaningless. It should be noticed that the 

asymptotic expansion result within the plateau region does not give the accurate value 

4 

of the energy. It differs from the iteration one in the order of e^s^, which provides the 
limitation of the accuracy of the asymptotic expansion. 

Recently it has been proved [5] that the /-iteration in one-dimensional problem is 
convergent if the trial function is properly chosen to have a finite perturbed potential 
w{x), satisfying the conditions w{x) > 0, w'{x) < and w{oo) = 0. There is no restriction 
to the magnitude ofw{x). For the double-well potential, Wey{x) defined in (3.10) satisfies 
these conditions when g' > 1. Our numerical results show that both /- and r-iteration 
are convergent for g > 1, although it is not an easy task to prove the convergence of the 
r-iteration. Furthermore, the two iteration series could also be applied in some region 
of (? < 1, where Wev{x) is still positive and finite but no more a monotonic function of 
X, and reasonable results are obtained. However, for very small g (e.g. g < 0.4), where 
w{x) becomes negative in some x-region, the /-iteration does not give reasonable results 
and the obtained En becomes unstable, while the r-iteration could still work, although 
the convergence becomes slower when g is smaller. The reason is the following: The 
iterative formula for the perturbed wave function / is expressed as sum of two terms. It 
could become negative if the term containing w{x) becomes negative. This would give 
a negative ip{x) in some x-region, while the solution for the ground state wave function 
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il>{x) should be always positive. For the r-iteration the perturbed wave function e"'^^^^ is 
always positive for any finite t(x), either positive or negative. This condition is fulfilled 
as long as is finite and ^^(oo) 0. Thus, the r-iteration gives less restrictions to 

the perturbed potential w{x). Therefore, it is of interests to further study the condition 
for the convergence of the r-iteration and to apply it to other physics problems where 
perturbation method could not be applied. 
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Appendix 



Appendix A. Asymptotic series expansion of t+ and 8+ 

Starting from the integral equation of t\_ of the first equation in (3.20) 

oo 

X 

an asymptotic series expansion of r+ and could be obtained. From (A.l) it is easy to 
obtain the following equation 

il 0^0' = [{u -S^) + l «)^]0^ (^.2) 

Using the definition of 0+ in (3.6), the above equation leads to 

gS',r[ = \rl- S[r'_, -(«-£+)-! (r;)^ {A.3) 

Now let us expand both t'^{x) and £+ in power series of as following 

cx) -1 00 -j 

and S^-j:-e^^,. (A4) 

1 y y 



Substituting (A. 4) into (A. 3), comparing terms proportional to g we obtain a series of 
equations for {Sm} and {em}'- 

g^ SqS'2 = -{u-ei) 

9-' S',Si = ls'i-S[S', + es-ls'i (A5) 

1 1 """^ 

^ ^ n=l 

The above equations are exactly the same as those obtained in Ref.[5], when taking 
S+ — —E and = —E^. Prom the first equation of (A. 5), considering u = ji^^i we 
have ei = I when setting x = 1. 

In the following we derive a compact expression of {S'^} and {e^}. Assuming, for 
m > 1, 

-1 2m-l 
^ i=0 
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where 

Considering S[ = from the last equation of (A. 5) we have 

1 2m- 1 

SoS'm+2 = E A(m)e'+'(/ + 4) 

1 m-1 2n-l 2(m-«)-l 
^ n=l i=0 j=0 

For the second summation on the right hand side of (A.8), defining i + j + 1 — I, we have 
Imin — 1 and Imax = 2m — 1, which leads to 

^ 2m- 1 



= E + 4) 



-1 2m— Im— 1 imox 

E E E 2A(n)A-i-i(m-n)e'+' + e^+i, (A.9) 



24m+2 

where 

^min = max(0, / — 2(m — n)) 
imax = min(2n - 1, / - 1). (^.10) 
To fix Cm+i we put X — 1, then ^ = 1 and S'q — 1 — — 0. This gives 

-1 2m— 1 1 2m— 1 m— 1 imax 

EAM(/ + 4) + ^^ S 2A(n)A-^-i(m-n). (All) 



Introducing Sq — — could be expressed as 

-j 1 2m— 1 

2m— 1 m— 1 imax 

+ E E E 2A(n)A-.-i(m - n)(r^ - 1)}. (^.12) 

/ — 1 Tl — 1 % — imin 

Applying the equality 

1+2 

_ 1 = _ 1)(^H2 + ^.+1 + . . . + 1) = _ 1) ^ (^_;L3) 

Z,=0 

changing the summation order of L and / in (A. 12) we finally reach the following expression 

of 'S'^+2- 
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-1 2m+l 2m— 1 2m— 1 m— 1 imax 

^'m+2 = ^^e E e"{E AM(/ + 4)+ E E E 2A(n)A-.-i(m-n)} 

Z,=0 (=Zl i=(2 "=1 i=imin 

-1 2m+l 

= ^4(;;^e' E/5a-+i)e^ (A14) 

^ L=0 

with 

Zi = max(0, L — 2), ^2 = niax(l, L — 2). 
For m = 1, it is easy to get /3o(l) = 1 and /3i(l) = 1. Finally we obtain , for m > 1, 

/3L(m + l) = /3°(m+l) + AA(m + l) (A15) 

where 

2 m- 1 

/3°(m + l)- E A(m)(/ + 4) (A16) 

;=max(0,L-2) 

2m— 1 m— 1 imax 

A/3i(m + l)= E EE 2A(n)A_,_i(m-n) 

i=max(l,L— 2) n=l i=imin 

Taking L = in (A.15) and (A.16), from (A.ll) we have 

Wi = ^i^/?o(m + l). (A17) 

Based on (A.15) and (A.16) a revised pyramid structure of (5i{m) could be constructed in 
a similar way as those in Appendix D of Ref.[l]. 

m — 1 



TO = 2 



m = 3 (A18) 

By using (A.15) and (A.16), we see that each row I3i{m + 1) can be obtained from the 
row I3i{m + 1) and APi{m + 1) while the later two could be obtained from Pi{m) above. 
For example, using /?i(l) = /5o(l) = 1, we have 























/3o°(2) 






A/33(2) 


A/32 (2) 


A/3i(2) 


A/3o(2) 










A (2) 


/3o(2) 








m 


/32°(3) 


/3?(3) 


/5o°(3) 


A/55(3) 


A/?4(3) 


A/33 (3) 


A/32 (3) 


A/3i(3) 


A/3o(3) 


/55(3) 


/34(3) 


/53(3) 


/32(3) 


/3i(3) 


/5o(3) 



/33°(2)=/3i(l)-(l+4) = 5, 
/32°(2) = /3i(l)-5 + /3o(l)-4 = 9, 
A/32 (2) = A/33 (2) = 
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etc. The values of the elements in the pyramid are 



m 








1 


1 








m 


= 1 








5 


9 


9 


9 










AA(2) 


























A (2) 






5 


9 


9 


9 






m 


= 2 






35 


89 


134 


170 


170 


170 








AA(3) 




2 


6 


8 


8 


8 


8 








A (3) 




37 


95 


142 


178 


178 


178 




m 


= 3 




333 


1093 


2087 


3155 


4045 


4757 


4757 


4757 






A A (4) 


20 


76 


148 


220 


256 


256 


256 


256 






A(4) 


353 


1169 


2235 


3375 


4301 


5013 


5013 


5013 


m 


= 4 



Correspondingly, the energies are 



9 178 89 5013 

^2=2^, ^3=2l^=2^, ^4 = ^, etc. (A20) 



Appendix B 

For the convenience of applying the coordinate system {S, a} defined in Eq.(2.10), 
the definition of some quantities introduced in ref.[2] is given in the following. For the 
new coordinate system 

{S, a) = {S, Q;i(q), Q;2(q), • • • , Q;jv-i(q)), {B.l) 
each point q in the iV- dimensional space will now be designated by 

{S,ai,a2, - ■ ■ ,oiN-i), (-B.2) 
instead of (gi, ^2, ^3, • " " > In)- The corresponding line element can be written as 

d Ci=s hsdS + ^ ttj hjdaf, {B.3) 



the gradient is given by 



Aid 1 9 
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The kinetic energy operator 

T= - 

2 



can be decomposed into two parts: 

T^Ts + T^, (B.e) 

with 

rp ^ ^_^('}^^) (B7) 

^ 2hsKdS^hsdS'' ^ ■ ' 

in which 

Af-l 

ha=I[ hj, (5.9) 

3=1 

and 

The volume element in the q-space is now 

d^q = hsKdSda {B.ll] 

with 

N-l 

da^ W day (-B.12) 

i=i 
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